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ABSTRACT 


Neubauer,  Albert  H.,  Jr.,  Ph.D.,  Purdue  University, 

June  1963.  The  Dynamic 8  of  an  Blaatle  Linkage.  Major  Pro- 
v  fessors :  A.  S.  Hall  Jr,  and  Raymond  Cohen. 

^^•As  the  operational  speed  and  accuracy  of  machinery  In¬ 
creases,  the  method  of  rigid  body  linkage  analysis  fails  to 
adequately  predict  the  dynamic  behavior  of  the  members.  As 
a  result  the  exact  position  of  a  member  at  any  given  time 
cannot  be  determined. 

•The  purporrc"  uf  this' •tnve'B t lgat  1  on  was  to  complete  an 
analytical  study  of  the  transverse  vibrational  character¬ 
istics  of  the  elastic  connecting  rod  la  an  in-line  slider 

.x,  ye 

-crank  mechanism^  A  non-linear  fourth  order  partial  differ¬ 
ential  equation  for  the  transverse  vibration  of  the  connect¬ 
ing  rod  has  been  derived.  This  equation  was  reduced  to  a 
linear  form  involving  coefficients  which  were  functions  of 
the  position  coordinate  on  the  connecting  rod- and  time. 

This  linear  equation  was  solved  by  a  finite  difference 
scheme  on  an  IBM  7090  digital  computer.  g&Z. 

In  addition,  a  further  simplifying  assumption  was 
applied  to  the  above  linear  equation  of  motion  which  reduced 
the  coefficients  to  functions  of  time  only.  This  simplified 
linear  equation  was  also  solved,  using  the  same  finite 


lx 


difference  scheme  as  above,  to  determine  the  validity  of 
this  assumption. 

The  simplified  linear  fourth  order  equation  of  motion  . 
was,  in  turn,  reduced  to  a  set  of  uncoupled  second  order 
equations  of  motion  by  making  use  of  the  technique  of 
normal  mode  expansions.  These  uncoupled  equations  were 
solved  on  the  IBM  7090  using  the  Runga-Kutta  method  of 
numerical  quadratures . . 

Correlation  between  the  various  solutions  was  achieved. 
The  solution  of  the  first  uncoupled  equation  of  motion  was 
found  to  be  a  periodic  approximation  of  the  finite  differ¬ 
ence  solutions. 

The  "two  dimensionless  parameters,  S  the  speed  parameter 
and  Z  the  physical  parameter,  that  naturally  arise  during 
the  derivation  of  the  uncoupled  equations  of  motion  can  be 
used  by  the  designer  as  an  Indication  of  the  vibrational 
activity  of  the  mechanism.  It  was  found  that  Z  must  always 
be  less  than  n*S  and  that  the  limiting  condition,  Z  =  »*S, 
is  to  be  avoided. 
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INTRODUCTION 

In  the  past  the  hulk  of  the  design  work  requiring  link¬ 
age  synthesis  has  been  accomplished  by  assuming  that  the 
members  or  links  were  rigid  bodies;  or  In  other  words  that 
no  elastic  deformation  takes  place  under  dynamic  loading 
conditions.  Thus  for  a  single  degree  of  freedom  system, 
such  as  a  slider-crank  mechanism  only  a  geometrical  relation¬ 
ship  between  the  input  and  the  output  is  required  and,  there¬ 
fore,  once  the  input  Is  specified  the  output  Is  completely 
determined.  If,  on  the  other  hand,  the  connecting  rod  Is 
assumed  to  be  elastically  deformed  because  of  tne  dynamic 
loads  inherent  in/body  motion,  then  the  relationship  men¬ 
tioned  above  becomes  Inadequate  to  describe  the  output  as  a 
function  of  the  input  and  other  relationships  are  required. 
The  output  will  however  execute  a  perturbative  motion  about 
the  predicted  rigid  body  position.  This  effect  would,  there¬ 
fore,  seem  to  present  a  problem  in  mechanism  design  if  these 
vibrations  were  large  enough  in  amplitude  to  cause  the  out¬ 
put  to  vary  beyond  some  arbitrarily  chosen  acceptable  limit. 
This  problem  is  especially  pointed  up  in  high  speed,  accurate 
machinery  where  dynamic  loads  are  high  and  members  have  low 
cross-sectional  area  to  length  ratios. 

It  will  be  shown  that,  in  general,  the  situation  reduces 
itself  to  a  complex  problem  in  beam  vibration  theory.  All 
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but  the  very  basic  vibrations  texts  treat  the  basic  theory 
of  the  free  transverse  and  axial  vibration  of  a  beam  and 

mathematical  solutions  are  presented  for  various  boundary 

1  # 

and  initial  conditions.  Some  more  advanced  texts  extend 
this  theory  to  consider  the  effect  of  axial  and  transverse 
loads  on  the  equation  of  motion  and  the  natural  frequencies 
of  the  system.  But  these  considerations,  although  important 
in  the  understanding  of  the  problem  defined  in  this  thesis, 
are  not  directly  applicable  to  the  synthesis  of  a  mechanism. 
Two  technical  papers  have  been  written  applying  vibration 
theory  to  mechanism  synthesis.  One,  dealing  with  the  longi¬ 
tudinal  vibration  of  the  coupler  of  a  crank-rocker  mechanism 

p 

was  written  by  Halit  Kosar  at  the  University  of  Istanbul. 
This  paper  derives  a  non-homogene ous  Mathleu  equation  to 
represent  the  differential  relationship  governing  the  varia¬ 
tion  in  output  crank  angle  as  a  function  of  the  input  crank 
angle.  Solutions  of  this  equation  of  motion  are  compared 
with  the  results  from  an  experimental  set-up.  The  other, 
written  ,by  Prof.  Dr.  Ing.  W.  Meyer  zur  Capellen^,  deals  with 
the  transverse  vibration  of  the  coupler  of  a  crank-rocker 
mechanism.  1  method  of  approach  quite  similar  to  that  used 
in  this  thesis  was  used  by  Prof.  Meyer  zur  Capellen  to  derive 
a  non-homogene ous  fourth  order  partial  differential  equation 
of  motion  with  non-constant  coefficients.  This  equation  of 


♦Numerical  superscripts  refer  to  entries  in  the  bibliography 
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motion  compares  with  equation  (20)  of  this  thesis,  and  re¬ 
presents  the  relationship  between  the  magnitude  of  the  trans 
verse  vibration  of  the  coupler,  the  position  along  the  coup¬ 
ler,  and  the  angular  Input  of  the  crank.  A  major  point  of 
divergence  between  this  investigation  and  that  of  Prof.  Meye 
zur  Capellen's  occurs  in  the  derivation  of  the  pin  reaction 
relationships.  The  later  of  the  two  authors  mentioned  above 
has  published  many  papers  on  subjects  closely  related  to 
problems  associated  with  elasticity  in  linkage  design  and 
is  probably  the  foremost  authority  on  this  subject. 

When  approaching  the  problem  of  vibration  in  the  con¬ 
necting  rod  of  a  slider-crank  mechanism  many  questions  nat¬ 
urally  arise  which,  in  the  course  of  an  investigation, 
should  be  answered.  In  this  case  one  of  the  first  questions 
is,  under  what  conditions  are  vibrations  in  the  connecting 
rod  important?  Also,  which  vibration  has  more  effect  on  the 
output,  transverse  or  longitudinal  vibrations?  If  the  in¬ 
vestigation  is  to  have  future  use  as  a  design  tool,  then  a 
relatively  simple  method  of  determining  the  effect  of  con¬ 
necting  rod  vibrations  on  the  output  is  required,  and  also 
the  effect  that  linkage  parameters  have  on  these  vibrations. 
To  have  confidence  in  the  validity  of  final  theoretical 
mathematical  results  some  correlation  between  theoretical 
and  experimental  data  should  be  achieved.  Finally,  the  re¬ 
sults  should  be  extended,  if  possible,  to  include  more  com¬ 
plex  data  such  as  input  functions  that  better  approximate 
actual  input  conditions. 
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Along  with  answering  some  of  the  above  questions,  the 
object  of  this  Investigation  Is  to: 

a)  Determine  the  equation  of  motion  for  the  transverse 
vibration  of  the  connecting  rod  of  a  slider-crank  mechanism 
In  Its  most  general  form. 

b)  Obtain  a  solution  to  the  equation  of  motion  by  re¬ 
stricting  the  problem  through  appropriate  simplifying  assump¬ 
tions. 

c)  Determine  the  effect  of  varying  linkage  parameters 


on  this  solution 
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ANALYSIS  OP  THE  SYSTEM 

Be script  Ion- of  the  Mechanism 

During  this  Investigation  only  an  ln-llne  slider-crank 
mechanism  will  be  considered.  See  Figure  1.  It  will  be 
assumed  that  the  connecting  rod  of  this  mechanism  is  a  beam 
that  obeys  the  fundamental  laws  of  elasticity.  For  the  pur¬ 
pose  of  simplification,  the  connecting  rod  cross-sectional 
area  will  be  constant  and  the  properties  of  the  rod  will  not 
vary  along  the  length. 

The  geometrical  equations  governing  the  motion  of  this 
type  of  slider-crank  mechanism  are  well  known,  when  all  of 
the  links  are  assumed  to  be  rigid  bodies.  See  Appendix  A. 

If  the  input  is  taken  to  be  the  angular  displacement,  e(t), 
of  the  crank  and  the  output  is  considered  to  be  the  position 
of  the  piston,  Xg,  then 

Xg  =  r  cos  eCt )  +  [La  -  r*  sin8  0(t)]^  (l) 

is  the  only  equation  required  to  determine  the  input-output 
relationship.  If,  on  the  other  hand,  the  connecting  rod  is 
assumed  to  have  elastic  properties  then  the  length,  L,  of 
this  rod  will  no  longer  be  constant  but  also  an  undetermined 
function  of  time.  The  above  equation  (1),  therefore,  is  no 
longer  adequate  to  completely  describe  the  system. 
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An  acceleration  analysis  of  the  rigid  body  mechanism 
can  easily  be  made  to  determine  the  components  of  accelera¬ 
tion  perpendicular  and  parallel  to  the  connecting  rod  for 
any  point  on  the  connecting  rod  (see  Appendix  A)  as  a  func¬ 
tion  of  time  and  displacement  along  the  rod.  At  the  point 
B,  or  at  x  =  L,  the  components  of  equation  (2a)  perpendicular 
and  parallel  to  the  connecting  rod  will  be  equal  to  the  de¬ 
rived  expressions,  equations  (3a)  and  (4a).  Using  this  meth¬ 
od  of  analysis,  equations  for  the  components  of  acceleration 
along  the  connecting  rod  and  perpendicular  to  the  connecting 
rod  can  be  obtained.  These  components  will  be  called  the 
rigid  body  components  of  acceleration.  However,  when  the 
assumption  of  rigid  body  mechanics  is  relaxed  acceleration 
components,  thus  derived,  do  not  represent  the  accelerations 
in  the  system  but  must  be  modified  to  Include  the  accelera¬ 
tions  due  to  vibration. 

If  a  general  point  C  on  the  connecting  rod  is  assumed 
to  vibrate  only  in  a  line  perpendicular  to  the  rigid  body 
position  O'  of  this  same  point  on  the  connecting  rod,  then 
from  Figure  2,  it  can  be  seen  that  the  position  of  this 
point  in.  terns  of  the  rigid  body  position  and  the  magnitude 
of  vibration  will  be, 

=  XQi  +  w  sin  <f> 

Yq  s  Yci  +  w  cos  <p 

The  first  derivative  with  respect  to  time  yields  the  velocity 


The  second  derivative  with  respect  to  time  yields  the  accel 
e ration 


Xc  =  Xci  +  if  sin  <p  +  2w  <p  cos  <p  +  w  <p  cos  <p  -  w  ipa  sin  <p 

•  •  *  •  #  t  ## 

Xc  =  Y0,  +  w  cos  <p  -  2 w  <p  sin  <p  -  w  #  sin  <p  -  w  cos  <p 

Transforming  from  the  (X,Y)  coordinate  system  to  the  (x,y) 
coordinate  with  the  orthogonal  transformation  equation* 


/  \ 

1 

(..  ) 

*c 

“  cos  <p  -  sin  <p- 

xc 

< 

= 

■< 

.. 

y'c 

_  8 In  <p  C08  <p_ 

Yc 

V  J 

1 

J 

yields  the  components  of  acceleration  parallel  and  perpen¬ 
dicular  to  the  connecting  rod 


xc  =  [Xci  cos  (p  -  Yci  sin  <p]  +  2w  ip  +  w  ip 
as  a(x,t)  +  2w(x,t)  <p(t)  +  w(x,t)  p(t) 


(4) 


yc  =  [Xqi  sin  <p  +  *Yc,cos  <p]  +  w  -  w  p3 

(5) 

=  0(x,t)  +  w(x,t)  -  w(xft)  ^(t) 

In  equation  (4)  a(x,t)  represents  the  rigid  body  accel¬ 
eration  of  the  general  point  C  in  the  x  direction;  2 tp  w 


*Care  must  be  exercised  when  using  this  coordinate  transforma¬ 
tion  because  the  transformation  matrix  is  a  function  of  time* 
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represents  the  Coriolis'  component  of  acceleration  and  w  ip 
represents  the  relative  tangential  acceleration  between 
points  C  and  C' .  In  equation  (5)  0(x,t)  represents  the  rigid 
body  acceleration  of  point  C  in  the  y  direction;  w(x,t)  rep¬ 
resents  the  vibrational  acceleration,  and  w  <pa  represents 
the  relative  normal  acceleration  between  points  C  and  O', 

See  Appendix  A  for  the  derivation  of  the  expressions  for 
ei(x,t),  equation  (3a),  and  /3(x,t),  equation  (4a). 

All  of  the  forces  acting  on  the  connecting  rod  and  pis¬ 
ton  are  shown  in  the  free  body  diagram  of  Figure  3«  The  pin 
reactions  at  the  left  end,  A,  and  the  normal  force  acting  on 
the  piston  are  unknown.  Also  unknown  is  the  effect,  of  the 
vibrational  acceleration  on  the  transverse  and  longitudinal 
inertia  components.  In  order  to  derive  expressions  for  these 
unknowns  as  well  as  the  eventual  derivation  of  the  equation 
of  motion  the  following  assumptions  were  made  : 

1 )  The  crank  and  piston  are  rigid  links  and  the  con¬ 
necting  rod,  link  3,  is  a  uniform  elastic  member. 

2)  All  of  the  physical  properties  of  the  mechanism, 
the  input  e(t),  and  the  force  on  the  piston,  P(t),  are  known 
quantities . 

3)  Friction  at  the  pins  A  and  B  and  on  the  slider  is 
assumed  zero. 

4)  The  vibration  of  the  connecting  rod  is  small  in  com¬ 
parison  to  the  length  of  the  connecting  rod. 

5)  The  vibration  of  the  connecting  rod  is  a  perturbation 
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/y,/3,w,w,w,etc. 


/ 

/ 


x 


Figure  4.  Sign  Convention  Diagram 


Q.  t 
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about  the  rigid  body  motion  of  link  3. 

6)  The  connecting  rod  will  vibrate  only  in  the  trans¬ 
verse  direction  and  the  transverse  acceleration  of  the  con¬ 
necting  rod  will  be  jj(x, t)  +  wCx,t).  See  equation  (5).  The 
normal  relative  acceleration,  w  ip2,  term  is  neglected  due  to 
assumption  4. 

7)  The  axial  acceleration  of  the  connecting  rod  is 
a(x,t)  only.  See  equation  (4).  The  Coriolis'  component  of 
acceleration,  2w  p,  and  the  tangential  relative  acceleration, 
w  ip,  terms  are  neglected  due  to  assumption  4.  4xlal  vibra¬ 
tion  is  also  neglected. 

Using  the  modified  expressions  for  the  longitudinal  and 
transverse  acceleration  of  link  3  the  longitudinal  and  trans¬ 
verse  inertia  terms  I  and  I  as  shown  in  Figure  3  can  now 

a  y 

be  defined  as : 

I_  =  JL  ms  «U,t)  d  i  (6) 

0 

Iv  =  JL  ma[£U,t)  +  wU,t)]  d  t.  (7) 

y  0 

Figure  4  shows  the  sign  conventions  that  will  be  used 
during  the  derivation  of  the  equation  of  motion.  In  general, 
the  positive  x  and  y  axes  define  positive  directions,  except 
that  positive  inertia  forces  are  directed  opposite  to  positive 
acceleration. 
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Derivation  of  the  Equation  of  Motion 
The  first  step  in  the  derivation  of  the  equation  of 
motion  is  to  determine  the  pin  reactions  Rx  and  Ry.  Making 
use  of  the  D'Alembert  principle  of  dynamic  equilibrium  and 
remembering  that  a  positive  inertia  force  is  in  the  direction 
of  negative  acceleration,  the  following  force  and  moment 
summations  can  be  made 

£  Fx  =  Rj.  “  Ix  -  N  sin  <p  +  (P-Ip)  cos  <p  =  0 ,  (8) 

£  Fy  =  Ry  -  Iy.  +  N  cos  «p  ♦  (P-Ip)  sin  (p  =  0.  (9) 

£  Mg  =  -L  Ry  +  jjj  m3  w(4,t)  a(4,t)  d  4  + 

pL 

JQ  ms(L-4)[0(4,t)  +  w(4,t)]  d  4  =  0.  (TO) 

From  (10)  the  reaction  Ry  can  be  determined  as 

R v  =  JL  -r  wU.t)  a(4,t)  d  i  +  JL  m9(1-  f)[^*(4,t)+w(4,t)]  d  4. 
*  0  0 

(11) 

Substituting  (11)  into  (9)  and  solving  for  N  gives 

Jl  Ite  }  C0(«,t)  +  Wt.tlJ  M  -  (P-I .)  sin  «],  (12) 

0  L  p 

Substituting  (12)  into  (8)  then  gives  the  remaining  pin 
reaction 
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\  =  ;L  ma[1-  xl Mift)  d  i  + 

0 

Iq  ms  £  tan  <p[0(4.t)  +  w(4,t)]  d  4  -  (P-Ip)  sec  (13) 

With  the  required  pin  reactions  R^  and  Ry  now  given  ae 
known  functions,  the  Euler-Bernoulli  equation 

-El  w"(x,t)  rr  L  M(x,t) 

can  be  employed  to  derive  the  equation  of  motion.  In  the 
above  equation,  the  negative  sign  appears  because  counter¬ 
clockwise  moments,  which  are  assumed  positive,  cause  negative 
curvature.  The  prime  notation  denotes  differentiation  with 
respect  to  x. 

Summing  moments  about  some  general  point  x,  see  Figure 
4,  and  substituting  this  summation  into  the  above  equation 
gives 

-  E  I  w"(x,t)  =  w(x,t)  -  Ry  x  + 

,fx  mg  (x-4)[/3(4,t)  +  w(4,t)]  d  4  - 

0 

JX  ms  [w(x,t)  -  w(4,t)]  a(4,t)  d  4  .  (14) 

0 

Substituting  for  Rx  and  Ry  the  expression  In  equations  (11) 
and  (13)  gives  the  complete  form  of  the  Buler-Bemoulll 


equation 
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-E  I  w"(x,t)  =  w(x,t)  [J^[1-  <]  o(4,t)  d  i  + 

m3  J  tan  <p[/3(4*t)  ♦  w(4.t)]  d  4  -  CP-Ip)  sec  v>]  - 
x[j^~wU»t)  ci(4,t)  d  4  +  J*^  m3(1--j*)[/3(4,t)+w(£,t)]d4] 

+  J*x  % (x-4)[0(4,t)  +  V(4,t)]  d  4  - 
0 

Jx  m3[w(x,t)  -  w(£,t)]  a(4»t )  d  4, 

0 

Rearranging  terms  for  convenience  yields  the  equation 
-E  I  w"(x,t)  =  -  w(x,t)[P-Ip]  sec  ip  - 

J*X  m3[w(x,t)  -  w(4,t)]  a( 4»t )  d  4  + 

J*  m,  Cw(i,t)  -  ifU.^1  taa  f  .  1  3  a({,tut 

+  Jx  m3(x-4)[/3(4,t)  +  w(  4»t)]  d  4  + 

0 

.  x(1.  ^)][^(4,t)  +  w(4,t)]  d  4.  05) 

In  an  attempt  to  eliminate  the  Integrals  In  the  above 
equation  C 1 5 )  the  standard  technique  of  using  Leibnitz*  rule 
for  the  differentiation  of  an  Integral  will  be  used.  There¬ 
fore  ,  differentiating  equation  (15)  with  respect  to  x  gives 
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-E  I  w'v(x,t)  =  -  w'(x,t)[P-Ip]  sec  i(>  - 

w' (x,t)  J* x  m3  a(l,t)  d  4  + 

0 

J*  %C»'(x,t)(t-  tout)  .  «U,t>  d  4  „ 

J*  ♦  w(i»t)]  d  4  ♦ 

0 

J*£  a3  .  (1.  £)][/3(4it)  +.w(4,t)]  d  4  . 

Differentiating  again  with  respect  to  x  gives 

-E  I  wIV(x,t)  =-w"(x,t )[P-Ip]  sec  <p  -  w"(x,t)  J*  n^a^t)  d  4 

+  w"(x,t)  m*[1-  tt]  a(4»t)  d  4  - 

w’ (x,t )[mja(x,t) ]  + 

w"(x,t)  rL  ms  *  [0(4, t)  ♦  w(4,t )]  d  4  + 

0 

mg[^(xtt)  +  w(x»t)].  (16) 

Without  looking  at  the  functions  which  represent  the  co¬ 
efficients  the  above  equation  (16)  can  be  written  in  the 
following  form 

wIV(x,t)  +  fx(x,t)  w"(x,t)  +  f,  (x,t)  w'(x,t)  +  A  V(x,t) 
w"(x,t)  J^[f3(4,t)w(4,t)  +  f*(4,t)w(4,t))  d  4  ♦  f*(x,t). 
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The  last  equation  can  be  classified  as  a  non-linear,  non 
-homogeneous,  fourth  order  partial  Integral-differential 
equation  with  non-constant  coefficients.  Equations  of  this 
type  do  not  lend  themselves  to  complete  mathematical  solu¬ 
tions  or  even  to  solutions  using  numerical  techniques. 

Within  the  frame  work  of  the  above  equation,  however, 
the  fourth  order  partial  differential  equation  governing  the 
free  transverse  vibration  of  a  beam 

wIV(x,t)  ♦  A  w(x,t)  =  0 

can  be  recognized.  The  additional  terms  arise  from  the 
forced  motion  of  the  beam  (connecting  rod)  through  its  con¬ 
strained  positions. 

Assumptions  for  LlnearlzatlOfl 
Since  the  derived  equation  of  motion  (16)  is  non-linear 
in  character  and  cannot  be  linearized  by  ordinary  mathematical 
techniques,  a  study  of  the  basic  components  of  the  Euler- 
Bemoulli  formulation  was  undertaken  to  determine  where  the 
non-linear  terms  arise.  By  following  the  various  components 
of  the  Euler-Bemoulli  formulation  (14)  through  the  process  of 
differentiation,  it  can  be  seen  that  all  of  the  non-linear 
terms  arise  from  the  derived  expression  for  R  ,  equation  ( 1 3 ) • 
R  ,  not  being  a  function  of  w  or  x,  drops  out  during  the  diff- 

y 

erentiation  with  respect  to  x.  A  study  of  Rx  was  therefore 
undertaken  to  determine  if  this  term  could  be  simplified. 

If  the  assumption  is  make  that  the  end  reactions  do  not 
depend  upon  the  vibration  then  all  of  the  non-linear  terms 
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are  neglected  and  the  pin  reactions  become  those  derived  If 
all  links  are  assumed  to  be  rigid.  By  eliminating  the  terms 
In  equation  (13)  that  depend  on  w  Its  time  derlvltlves  a 
relatively  complex  function  of  time  Is  obtained 

Rx=  JL  “»  d  t  +  + 

X  0 

fL  m3  £  tan  <p  [(}(  4»t )  ]  d  4  **  (P-Ip)  sec  (p .  ( 1 3' ) 

0  L 

This  function  (13*)  can  be  substituted  Into  equation  (14) 
and  after  differentiating  twice  with  respect  to  x  a  simpli¬ 
fied  form  of  equation  (16) 

-B  I  wIV(x,t)  =  -  w''(x,t)[P-Ip]  sec  <p  - 


wM(x,t)  fx  m3  «(€*t)  d  4  + 

0 

w"(x,t)  rL  ms  a(4,t)  d  4  -  w'(x,t)  m*  a(x,t)  + 

0 

w"(x,t)  JL  -3  *  /}(  4,t)  d  4  + 

0  L 


maCjsfx.t)  +  w(x,t)] 


(16') 


is  obtained.  This  equation  ( 1 6 * )  is  a  non-homogene ous, 
linear,  fourth  order  partial  differential  equation  with  non 
-constant  coefficients  that  can  be  solved  on  a  high  speed 
digital  computer  by  using  numerical  techniques. 

However,  in  the  hope  of  further  reducing  the  complexity 
of  equation  ( 1 6 * )  the  rigid  body  force  analysis  of  Appendix 
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B  was  undertaken.  The  result  of  this  analysis  is  the  assump¬ 
tion  that 

Rx  =  -  K  cos  0 ( t ) .  (17) 

Where  K  Is  a  positive  constant  that  can  be  approximated  by 

K  =  (M,  +  M4)  r  ufi  .  (18) 

It  must  be  noted  that,  In  order  to  make  this  simplifying 
assumption,  the  scope  of  the  Investigation  was  restricted 
to  mechanisms  with  small  r/L  and  Ma/M*  ratios  and  to  input 
functions  that  are  linear  in  time,  0=ut.  If  these  restrict¬ 
ions  are  unsuitable  then  equation  (16')  must  be  solved  directly. 

A  further  simplifying  assumption  arises  from  the  in¬ 
spection  of  the  last  term  in  equation  (14).  In  light  of  as¬ 
sumption  4  on  page  9,  it  can  be  noted  that  the  term 

w(x,t)  -  w(4,t) 

is  proportional  to  the  slope,  w' .  Therefore,  in  light  of 
the  required  assumption  made  during  the  derivation  of  the 
Euler-Bernoulli  equation  that  the  slope  must  be  small,  the 
term  written  abovfe  can  be  neglected.  Assuming  that  this  term 
equals  zero  should  not  introduce  into  the  problem  an  error 
of  magnitude  greater  than  that  introduced  previously  by 
assuming  that  the  slope  is  zero.  This  assumption,  therefore, 
keeps  the  right  and  left  hand  sides  of  equation  ( 1 4 )  within 
the  same  order  of  magnitude. 
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Derivation  of  Simplified  Equations  of  Motion 
If  the  function  representing  Rx  in  equation  (17)  is 
substituted  Into  the  Euler-Bernoulli  equation(l4)  then  the 
equation 

-E  I  w"(x,t)  =  -  w(x,t)  K  cos  0  +  x  Ry.  + 

JX  ms  (x-4)[£(  4,t)  +  w(4,t)]  d  i  - 
0 

Jx  m3[w(x,t)  -  w(4,t)]  ci(C,t)  d  4  (19) 

0 

results.  Employing  the  same  technique,  equation  (19)  is 
differentiated  twice  with  respect  to  x 

-E  I  w"'(x,t)  =  -  w'(x,t)  K  cos  0  +  Ry  + 

JX  ms  C/3(  4»t )  +  w(4,t)]  d  i  -  w'(x,t)  Jx  ma  a(4,t)  d 
0  0 

-E  I  wIV(x,t)  =  -  wM(x,t)  K  cos  0  +  ma[0(x,t)  +  w(x,t)]  - 

w"(x,t)  fx  ms  a(4,t)  d  (  -  w'(x,t)Cm3  a(x,t)].  (20) 

0 

Equation  (20)  is  the  linear  form  of  the  equation  of  motion. 

If  the  assumption,  that  the  last  term  in  equation  (19) 
can  be  neglected,  is  now  included,  equatlon(19)  becomes 

-E  I  w”(x,t)  =  -  w(x,t)  K  cos  0  +  x  Ry  + 

Jx  ms (x»4)[£( 4,t)  +  w(4,t)]  d  4.  (21) 

0 
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Again  employing  the  same  technique,  equation  (21)  Is  differ¬ 
entiated  twice  with  respect  to  the  displacement  x  and  the 
equation 

-E  I  wIV(x,t)  =  -  w"(x,t)  K  cos  8  +  ms  [/3(x,t)-*-w(x,t)  ]  (22) 

results.  This  equatlon(22)  Is  the  most  simplified  form  of 
the  equation  of  motion.  With  the  exception  of  the  middle 
term,  equation  (22)  Is  the  equation  of  motion  for  the  trans¬ 
verse  vibration  of  a  beam  with  an  external  forcing  function. 
The  middle  term  In  equation  (22)  represents  the  effect  of  a 
time  varying  axial  end  load  on  the  transverse  vibration  of  a 
beam. 

To  complete  the  statement  of  the  mathematical  problem 
four  boundary  said  two  Initial  conditions  must  be  specified. 
The  boundary  conditions  are  the  same  as  those  for  a  beam 
with  pinned  ends.  In  terms  of  the  symbols  used  these  axe 

w(0,t)  =  0,  w"(0,t)  =  0, 

w(L,t )  s  0,  w"(L,t)  =  0. 

These  conditions,  stated  In  words,  are:  The  deflection  of 
the  beam  ends  Is  zero  and  the  moment  (second  derivative)  at 
the.se  ends  Is  also  zero.  If  the  assumption  Is  made  that  the 
connecting  rod  Is  started  from  0=0°  from  rest,  the  Initial 
conditions  are 


w(x,0 )  =  0, 


w(x,0)  =  0 
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These  conditions,  stated  in  words,  are:  The  connecting  rod 
is  initially  a  straight  beam  and  the  motion  is  started  from 
a  rest  position,  no  initial  velocity. 


Reduction  of  the  Equation  of  Motion  to  a  Set  of  Uncoupled 

Equations 

A  standard  mathematical  technique  for  obtaining  the  solu¬ 
tion  to  a  partial  differential  equation  is  the  method  of 
separation  of  variables.  To  use  this  technique  the  assump¬ 
tion  is  made  that  the  solution  is  of  the  fom 


w(x,t)  =  £ 
n=1 


qn(t)  <pn(x) 


(23) 


where  qn(t)  is  only  a  function  of  time  and  x )  is  only  a 
function  of  displacement.  The  assumed  solution,  equation 
(23),  is  then  substituted  into  the  partial  differential 
equation  and  the  equations  governing  the  functions  of  dis¬ 
placement,  9n(*)>  and  the  functions  of  time  qn(t),  are  deter¬ 
mined.  To  mathematically  determine  these  functions  it  is 
-  . 

necessary  to  prove  that,  for  a  vibrating,  system  of  this  type, 
the  modes  of  vibration,  ^n(x),  are  orthogonal,1  that  i3 


J*q  *n(x)  «>m(x)  =  0 


m  ^  n 


(24) 


It  has  been  shown  in  the  literature  that,  for  the  free 
transverse  vibration  of  a  beam  with  pinned  ends,  the  modes 
of  vibration  are1 


<pn(x)  =  sin 


(25) 
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The  trigonometric  functions  are  known  to  be  a  complete  set 
of  orthogonal  functions;  therefore,  the  conditions  stated  in 
(24)  are  fulfilled.  It  can  be  noted  also  that  the  <pn(x) 
functions  satisfy  the  boundary  conditions  as  stated  pre¬ 
viously. 

For  the  solution  of  equation  (22)  the  assumption  will, 
therefore,  be  made  that  the  solution  is  of  the  form 

w(x,t)  =  L  q  (t)  sin  (26) 

n=1  n 

and  the  time  functions,  q  (t),  will  be  determined. 

n 

Substituting  equation  (26)  and  its  required  derivatives 
with  respect  to  both  time  and  displacement  into  the  equation 
of  motion  (22)  gives 

B  I  2  q^t)  (fl5>)4  ♦  “3  £(x,t)  + 

n=1  ^ 

m,  S  q  (t)  sin  + 

n=1  n  L 

K  cos  6  qn(t)  sin  ^  =  0  .  (27) 

Making  use  of  the  orthogonality  condition  by  multiplying 
the  above  equation  by  sin  mrrx/L  and  integrating  over  the 
length  of  the  connecting  rod,  yields  the  relationship,  a 
second  order  ordinary  differential  equation,  that  the  qn(t) 
must  satisfy.  Performing  this  operation  yields 
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«  j,  *n(t>  ttT  .m  up- .in  dx* 

^  m3  j?(x,t)  sin  dx  +  Ba  J*q  E  '^(t)  sln®J*  sin  dx  + 

K  cos  9  jJ  l  ,n(t)  [a§J  .in  ®|*  sin  «P  dx  .  (28) 
n=1 

Interchanging  the  sum  and  integral  signs  and  noting  that, 
from  orthogonality 

J*  sin  *£*  sin  »}*  dx  =  0  m  ji  n 

and 

•1*0  slaa  “P  dx  =  2 

yields  (all  terms  except  when  n=m  are  dropped  from  (28) 


-ma  J“q  j3(x,t)  sin  B£i  dx 
or  for  convenience 

•i  (t)  ♦  [«<»)♦  „  iat)‘  K  sm  t]  ,  (t)  = 

n  LL*m3  L»m3  J  n 

-  £  Jo  2(x»t)  sln  ®P  dx  (29) 
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Since  0(x,t)  is  a  known  function,  equation  (4a'),  the 
integral  on  the  right  hand  side  of  equation  (29)  can  be 
evaluated  as  a  function  of  time  only. 

Bquations  (29)  are  known  as  the  uncoupled  equations 
of  motion  and,  in  this  case,  are  in  the  form  of  non- 
homogeneous  Mathieu  equations.  The  natural  frequencies 
for  the  various  modes  of  vibration  are  represented  by 
the  term 

w  a  _  [ LoaJt  SI  +  (ng)a  K  cos  8~| 
n  L  L*m,  L9m3  J  * 

It  can  be  noted  that  the  natural  frequencies  of  this 
system  vary  with  time.  This  is  to  be  expected  since  the 
axial  loads  are  time  varying  functions. 

To  further  verify  the  validity  of  equations  (29)  as 
the  uncoupled  equations  of  motion,  the  Lagrangian  approach 
is  employed  to  derive  these  equations. 

Making  use  of  the  same  assumptions  that  were  required  to 
arrive  at  the  equation  of  motion  ( 22 )  the  kinetic  and 
potential  energies  can  be  formulated  as  follows : 

1  =  7  CV  +  7C*  ]  d* 

where  x  and  y  are  defined  by  equation  (2)  since  during  an 
c  c 

orthogonal  transformation  the  length  of  a  veotor  remains 
constant.  Therefore, 


♦  y  9  =  X  * 

*  C  0 


+  Y  8 
*  Lc  • 
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The  transformed  set  of  equations  for  the  velocity  of  the 
general  point  0  are  not  used  because  the  transformation 
matrix  Is  a  function  of  time  and  later,  when  substituting 
this  kinetic  energy  term  into  Lagrange  s'  Equation,  dif¬ 
ferentiation  with  respect  to  time  la  required  and  a  term 
will  be  lost.  The  kinetic  energy,  therefore,  is 


T  = 


2~  So  C<Xc,  +  w  sin  <p  +  w  ip  cos  <p)*  + 


/  *  *  *  \ 

(Y  ,  +  w  cos  <p  -  w  <p  sin  <p)s]  dx 
c 

=  ~  Jo  CXc,»  +  Ye,a  +  2w  (Xc,  sin  <p  +  Yc,  cos  <p)  + 


2w  Xc  t  cos  <p  -  Yct  sin  <p)  +  it®  ♦  w®  <pa  ]  dx  (30) 

By  inspecting  two  of  the  terms  in  the  above  equation  it 
can  be  seen  that  the  coefficient  of  2w  is  the  component 
of  the  rigid  body  velocity  of  point  O'  in  the  y  direction 
and  the  coefficient  of  2w  <p  is  the  component  of  the  rigid 
body  velocity  of  point  O'  in  the  x  direction. 

If  the  form  of  the  solution  is  again  assumed  to  be  the 
infinite  sum  of  the  product  of  a  function  of  time  and  a 
function  of  displacement  then  the  expression  (26) 
can  be  substituted  into  the  kinetic  energy  equation  (30). 

In  this  situation  the  <ln(t)  are  thought  of  as  an  infinite 
set  of  generalized  coordinates  that  completely  describe 
the  motion  of  the  system. 
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Substituting  this  assumed  form  of  the  solution  and  its 
appropriate  derivatives  into  equation  (30)  yields  the 
final  form  of  the  kinetic  energy  term. 

T  =  t3Jo  IV5  ♦  V’  + 

cm  # 

2  E^  qn(t)  sin  (xc'  3*n  ^ o'  008  ^  + 

cm 

2  E  qn(t)  sin  ap 

C<p(icc,  cos  <p  -  Yc,  sin  <p)]  ♦ 
n=1 

E  q.n8(t)  sin*  ♦  V>a  E  qn8(t)  sin8  ap  )  dx. 
n=1  n=1 


The  potential  energy  can  be  written  as  the  sum  of 
the  potential  energy  due  to  bending  and  the  potential 
energy  due  to  the  axial  load.1 

V  =  *§  Jo  0"(x,t)]8  dx  +  J*  [w*  (x,t)]8  dx. 

Substituting  the  first  and  second  derivatives  with  respect 
to  x  of  the  assumed  form  of  the  solution  into  this  equation 
for  the  potential  energy  yields 

v  ■  S  Jo  j,  v<*> 1 OT  •«-  “P  »*  * 

*  vu>  W  »«*  “P**. 

n=  1 
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Interchanging  the  summation  and  Integral  signs  and  evaluating 
the  integrals  reduces  this  equation  to 

V  =  3ilL_  g  n*  q  *(t)  +  H.  &  g98..&.  £  n»  q  8(t). 

4L3  n=1  n  4L  n=1  n 

Lagranges'  Equations  for  this  situation  take  the 
general  form 

ft  [sf^TTT]  -  +  *j|¥tT  =  0 

The  first  term  becomes 

dtL^rt)]  =  It  [“3  lo^c'  8ln  V  *  V  cos  sin  ap  dx  + 
4n(t)  Jo  m3  sin8  dx] 

sms  Jo  Hxc>  sln^)  +  Xci  <p  cos  sP  ♦  7C»  cos  <p  - 
Yci  ip  sin  <p)  sin  ^P  +  qn(t)  sin8  apj 

dx  0 

By  evaluating  the  integral  wherever  possible  and  by  re¬ 
placing  those  terms  that  represent  the  transverse  component 
of  the  rigid  body  acceleration  with  /3(x,t)  the  above  be¬ 
comes 

It  [  aft'l ) ]  =  m3  Jo  ^(x»t}  8ln  “f*  dx  *  “F*  <*n(t)  + 

n 

ms  Jq  4>^ci  cos  <p  -  Yc,  sin  <p)  sin  ■P  dx  • 
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The  second  term  becomes 


=  m3  Jo  «P(XC,  cos  (p  -  Yc,  sin  <p)  sin  *£*  dx  + 
m3Lfra 


2  Hn 


The  third  term  becomes 


_M_  =  Slim  ,  (t)  ♦  a8!*  £.  sag  a  ,  (t) . 

2l»  n  SL  “ 


Combining  these  terms  yields  the  uncoupled  equations  of 
motion 

~jT  *4n(t)  +  m3  Jq  |8(x,t)  sin  flj1  dx  + 

m3  Jq  £(Xc,  cos  p  -  Yc,  sin  <p)  sin  dx  - 
m3  Jq  p(Xci  cos  (p  -  Yci  sin  <p )  sin  apax- 


m3L  <pa 
2 


Rearranging  tenns  for  comparison  purposes  and  canceling 
the  third  and  fourth  terms  yields 

q.(t)  *  [a^LSI  t  Cli  E.M»  8  .  i,l  .  (t)  =  . 

L4  m3  L*  m3 


i  Jo  ®in  *5*  dx 


(33) 
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Equation  (33)  agrees  vlth  equation  (29)  term  for  term  except 
for  the  third  term  In  the  coefficient  of  qn(t).  This  term 
does  not  appear  In  equation  (29)  because  of  the  assumption 
that  the  normal  relative  acceleration  term  In  the  trans¬ 
verse  component  of  acceleration  can  be  neglected.  See 
assumption  6  on  page  1 1 . 
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SOLUTIONS  OF  THE  EQUATIONS  OP  MOTION 


As  derived  in  the  preceding  section  the  uncoupled  equa¬ 
tions  of  motion  are 


m 


£L 


qn(t)  =  - 


^  0(x,t)  sin  dx.  (29) 

Substituting  the  derived  expression  (4a1)  for  0*(x,t)  into 
equation  (29)  and  evaluating  the  integral  on  the  right  hand 
side  reduces  this  equation  to 


(t)  +  [ 


n*  it*  El 

l>mT 


2g*(t)  f  . 

“3? —  L1-008  n*\ 


2L  g*  (t ) 

m 


-cos  ntr 


J. 


(34) 


Where  gg(t)  and  g*(t)  are  defined  in  equations  (4a*)  and 
(4an)  of  Appendix  A;  K  is  defined  in  equation  (18),  and 
n  =  1,2,3,* • * . 

When  n  is  an  odd  integer,  n  =  1,3,5,***,  equation  (34) 


reduces  to 
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v*>  *  w-  -  aijVieMJl3  «.«*>  =  - 

55-  [s3(t)  +  \  g*  C-b  )3 


(35) 


where  the  right  hand  side  In  proportional  to  the  transverse 
component  of  acceleration  of  the  oenter  of  the  connecting 
rod. 

When  n  Is  an  even  Integer,  n  =  2,4,6,*- •,  equation  (34) 
reduces  to 


2L  g^  ( t ) 
nw 


(36) 


where  the  right  hand  side  Is  proportional  to  the  relative 
tangential  acceleration;  Substituting  the  expressions  for 
ga  (t),  g4(t)  and  K  into  equations  (35)  and  (36)  yield  the 

a 

complete  form  of  the  uncoupled  equations  of  motion 


>ir*BI 


ra.4j..  -  . 

L  wm,  + 


n*  ir 8  ( 1  + 


n=1 ,3,5*  * 
n=2 ,4, 6-  * 


where  Qq  =  ~  r  oj^cos  8  sin  <p  +  sin  6  cos  <p)  - 


L 


Slfl  S  .  L IWd  cos»8  sin 


cos  <p  L*  cos3<p 


*)] 


and  Q„  = 


Zk 

nir 


flfcd  Blfl  fl  .  lifcli. 

L  L  cos  ^ 


aaalfl  aJLa 

COS3^) 


*]. 


(37) 

(38) 


Since  8,  the  angular  input,  is  a  more  desirable  variable 
than  time  and  since  there  exists  a  simple  relationship  between 
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6  and  t,  0  s  to  t,  a  change  of  variable  from  t  to  0  In  equa¬ 
tions  (37)  and  (38)  yields  after  rearranging  terms  on  the 
right  hand  side  and  remembering  that 


d*qn  da qn 

dt3  d08 


Ve)  +  [££$*  §  oob  e]  «„<*) 

where  ^  ^2(cos  0  sin  <p  +  sin  0  cos  <p)  - 


ala  fl  +  £  Qoa*0  sin 

COS  (p  L  COS3  (p  J 

and  Q'  =  2L  fsin  6  _  £  fiflfl8  fl  ftlfl  tfll 

©  hit  Loos  <p  L  cos3  J  * 


(37’) 


(38*) 


The  dot  notation  now  denotes  differentiation  with  re¬ 
spect  to  0. 

Inspection  of  the  terns  in  equations  (37')  and  (38') 
shows  that  two  dimensionless  ratios  naturally  arise  that  in¬ 
clude  all  of  the  parameters  of  the  system.  These  dimension¬ 
less  equations  cure 


I^SgO)3 

(39) 

(40) 

S  can  be  thought  of  as  the  dimensionless  parameter  primarily 
describing  the  input  angular  velocity  and  Z  can  be  thought  of 
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as  the  parameter  that  describes  the  mechanism. 

The  right  hand  sides  of  equations  (37')  and  (38')  have 
not  been  non-dimenslonalized  but  depend  upon  the  parameter 
r,  the  length  of  the  driving  crank,  and  the  dlmftislonless 
ratio  r/L.  These  equations,  however,  could  be  non- 
dimenslonalized  by  redefining  the  vibrational  amplitude  as 
5n(0)  =  qn(e)/r. 

The  problem  would  then  reduce  itself  to  a  study  of  the 
effects  of  varying  the  three  dimensionless  parameters 
r/L,  Z  and  S.  Since  Z  depends  upon  the  r/L  ratio,  the  effects 
of  varying  r/L,  Ma/M*  and  3  on  the  vibrational  activity  of 
the  mechanism  would  be  sufficient. 
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and 


f  (e)  =  &  -  *  g-oaa  fl.ala  b~|  n  =  2  A  6. 

n  ntr  Lcos  <p  L  cos*  <p  J  *  * 


Equations  (41  )  in  conjunction  with  initial  conditions  are 
solved  simultaneously  to  obtain  the  deBired  result,  qn(t). 
Since  it  has  been  assumed  that  the  mechanism  will  be  started 
from  rest,  and  since  qfl(t)  is  the  function  that  describes  the 
behavior  of  time,  the  initial  oonditlons  for  equations  (57') 
and  (38' )  are 

Q.n(o)  =  4n(o)  =  0 


or  for  equations  (41 )  are 

V(0)  =  qn,(0)  =  °- 


Bffect  of  Damning 

The  above  initial  conditions  may  not  be  desirable  in 
some  future  investigation  of  this  problem.  If,  however, 
damping  can  be  introduced  Into  the  equations  of  motion  only 
the  steady  state  solution  will  remain  after  the  transients 
caused  by  free  vibration  have  been  damped  out.  The  uncou¬ 
pled  equations  of  motion  are  the  only  form  of  the  equations 
of  motion  thus  far  derived  that  yield  a  relatively  simple 
opportunity  to  study  the  effects  of  damping.  These  equa¬ 
tions  of  motion  are  in  the  well  known  uncoupled  form 


qn(t)  ♦  q^Ct )  =  f(t) 
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which  is  also  the  equation  of  motion  for  a  simple  spring 
-mass  system  under  the  influence  of  a  forcing  function.  To 
add  viscous  damping  to  this  equation  it  is  necessary  to  in¬ 
clude  a  term  that  is  proportional  to  the  damping  coefficient 
multiplied  by  the  velocity  (first  derivative).  Thus  the 
well  known  equation  for  the  motion  of  a  simple  spring-mass 
-dashpot  system  under  the  influence  of  a  forcing  function  is 
obtained. 

'^(t)  +  2tj  wn4n(t)  +  wa  qn(t)  = 

where  n  is  the  damping  ratio  (actual  damping  coefficient/ 
critical  damping  coefficient).  Since,  in  this  case,  the 
value  of  t)  is  completely  unknown,  various  different  values 
of  n  must  be  tried  until  the  steady  state  response  is  de¬ 
tected  in  the  solution. 

To  observe  the  effect  of  damping  on  this  system,  there¬ 
fore,  the  following  system  of  first  order  differential  equa¬ 
tions  was  programed  into  the  digital  computer, 

kr,  (0)  =  -  2ij  [n*ir*S  +  n*ir®2  cos  0]^q  (0)  - 

[n*ir4S  +  nairaZ  cos  0]  ( 0 )  +  fn(0)  (42) 

where  the  fn(0)  are  described  above.  The  Runga-Kutta  method 
of  numerical  quadratures  was  also  used  to  solve  the  set  of 
equations  (42). 
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Partial  Differential  Equations 
In  proceeding  from  one  form  of  the  equation  of  motion 
for  the  transverse  vibration  of  the  connecting  rod  of  a 
slider-crank  mechanism  to  another  form,  simplifying  assump¬ 
tions  were  required.  In  order  to  determine  the  validity 
of  some  of  these  assumptions  the  fourth  order  partial  dif¬ 
ferential  equations  of  motion  (20)  and  (22)  were  solved. 

The  high  speed  digital  computer  IBM  7090  was  again  used  for 
the  computation  and  a  FORTRAN  program  was  compiled  using  an 
eleven  point  finite-difference  approach. 

By  collecting  terms  and  substituting  the  representative 
equations  for  K,  expression  (18),  a(x,t),  expression  (3a1) 
and  0(x,t),  expression  (4a' ),  equation  (20)  can  be  rewritten 
as 

M 

w*V(x,t)  +  w(x,t)  -  [(1  +  j^*)  r  L  u?  cos  0  - 


x  r  u?(cos  0  cos  <p  -  sin  0  sin  <p)  -  f  W  ]• 

^w''(x,t)J  +  £r  u?  (cos  0  cos  <p  -  sin  0  sin  tp)  + 


c,og!. 

JLI 

l  L  ) 

cos* 

<pj 

r  of  (cos  0  sin  <p  +  sin  0  cos  <p)  - 

zrJLJif  3 la  fl  .  ffgSa,  a-gteLflh 

1  L  COS  <p  V  L  J  COS3  (0  J 


To  present  the  above  equation  in  a  more  desirable  form  the 
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change  of  variable  from  t  to  0  and  from  x  to  z  Is  required, 
where  0  =  u>  t  and,  therefore, 

i*w£x.t)  _  wa  g  JfUiflJ 

dt*  d  0a 


and  where 


L  VS 


so  that 

If  =  rW  ft*  =  rhs 

S  Is  the  dlmenslonles8  parameter  defined  In  equation  (39)* 
Performing  this  operation  yields  the  equation 

w*V(x,0)  +  w(x,0)  -  cos  0  -  (£)  (cos  0  cos  <p  -  sin0sirup) 

-  (£)a  SSfiifl.]  w"(x, 0)  +  [-*-  (cos  0  cos  (p  -  sin  0  sin  <p) 

2  h  cosa<pJ  lVSL 

+  x  (£f  CMa  A,!  v'  ( x,  0 )  = 
h  cosa  <pJ 

r(cos  0  sin  <p  +  sin  0  cos  <p)  - 


x  vs  {r  aliL-i  _  xl  aaal  fl  sin  a 3 

COS  <P  L  CC  !J3 


(43) 


In  this  case,  the  prime  notation  now  denotes  differentiation 
with  respect  to  x  and  the  dot  notation  denotes  differentia¬ 
tion  with  respect  to  0*  It  can  be  seen  that  this  differential 
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equation  of  motion  also  depends  upon  the  same  parameters, 

S,  Z,  r  and  r/L,  as  the  uncoupled  equations  of  motion  pre¬ 
viously  derived.  The  dependence  of  this  equation  on  the 
parameter  r  can  be  eliminated  by  redefining  the  vibrational 
amplitude  as 

w(x,0)  =  w(x,0)/r 

then  the  three  dimensionless  parameters  S,  Z  and  r/L  can  be 
varied  to  determine  their  effect  on  the  vibration  of  the 
connecting  rod. 

To  complete  the  statement  of  the  problem,  the  boundary 

conditions  for  the  solution  of  equation  (43)  written  in 

terms  of  x  and  0  are 

w(O,0)  =  0,  w( — 1-  ,  0)  =  0, 

Vs 

w*1  (0,0)  =  0,  vM( — —  ,  0)  =  0, 

Vs 

and  the  Initial  conditions  are 

w(x,0)  =  0,  v(x,0)  =  0. 

To  simplify  the  notation  of  equation  (43)  let 


f(x,0)  s  ~  cos  0  +  -2-  (£)  (cos  0  cos  -  sin  0  sin  <P)  + 

/S  Vs  L 

2  *1/  cos*  <PJ 

f'(x.fl)  =  —  (cos  8  co 8  <P  -  sin  0  sin  <P)  +  x  (£) 

%/SL 


and 


g(x,0)  =  [~r(cos  0  sin  o  +  sin  0  cos  <o)  - 

x  vs  {r  slnJL  -  4  I. 

cos  0  L  cos"  0  *  3 
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Therefore  equation  (42)  can  now  be  written  In  the  condensed 
form 

wIV(x,0)  +  w(x,0)  +  |~f(x,0)  w'  (x,0)  j  =  g(x,0)  (44) 

By  arranging  the  grid  and  using  the  difference  scheme 
as  shown  in  Figure  5  the  following  difference  equation  was 
derived^  to  represent  equation  (44) 

*3+1,11+1  +*3-i  ,k+1  "(a  *  1a  fi»l*3,k+1  = 

2(1  -jj-fjt  -  (Ax)1  f3,k)  *3,k  + 

(-a  +  (A  x)a  f3jk  +  iV*s  fj>k)  + 

(-a  +  (Ail*  f3ik  »3.,,k  + 

*3+a,k  *  *3-2, k  *  (2  +  W1  *3,k-1 

*3+1, k-i  -  *3-1, k-i  -  «3, k  (45) 

In  this  difference  equation  3  represents  the  increment 
of  displacement,  Lx,  along  the  connecting  rod  and  varies 
from  0  to  n+1 ,  the  total  number  of  displacement  intervals; 
and  k  represents  the  increment  in  angular  input,  L  0,  that 
varies  from  time  equals  0  to  any  arbitrary  stopping  time. 
The  boundary  conditions  reduce  to 


wO,k  -  °» 

wn+1 ,k  =  °» 

w1,k  =  "w-1,k, 

wn,k  =  "wn+2 
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Figure  5.  Grid  for  Finite  Difference  Solution 
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and  the  initial  conditions  become 

*3,0  =  °>  *3.1  *  H3.-f 

To  program  this  equation  into  the  computer  the  recursion 
matrix  form  of  equation  (45)  is  useful.  This  can  be  written 
as 

A  vrk+1  =  Wk  -  A  VT^"1  -  Gk  (46) 


where  the  matrix  A  is  the  following  (n  x  n)  symmetric,  trl 
-diagonal  matrix  with  constant  elements  that  depend  upon 
the  displacement  and  time  increments  only 
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Matrix  Is  the  (n  x  n)  matrix  of  coefficients 

Cyf-1  )  C6?+c?5  1  0 

(fiS-cf)  yS  (df*JF)  i 

1  yf  (6?+cJ) 

B  o  1  (fif-cj)  y£ 

e  •  e  e 

•  •  •  • 

k  ,  a&  B4  ^  f  ) 

where  y  =  (2 - -  2(Ax)aIJ,k;* 

(A  0)a 

6*  =  (-2  +(Ax)af3>k), 

(Ax)3 

and  <3  =  2  f  J»k* 

V 

The  vector  G  represents  the  non-homogene ous  or  right  hand 
side  of  the  equation  of  motion  (43) 

Gk  =  (A  x)4[g1  jk,  g2>k,  g3>k,  8n>kJ. 

It  can  be  seen  from  the  matrix  equation  (46)  that  the 
solution  of  the  difference  equation  at  any  given  Increment, 
w^*1 ,  of  time  depends  only  upon  the  values  of  the  function 
at  the  two  preceding  time  Increments,  w^  and  w^”1 ,  and, 
therefore,  In  solving  this  equation  the  behavior  of  the 
function  at  successively  Increasing  time  Increments  Is 
determined  In  a  systematic  manner. 

If  the  same  procedure  Is  used  for  the  finite-difference 
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solution  of  the  equation  of  motion  (22)  as  was  used  on 
equation  (20),  then  the  function  notation  used  in  equation 
(44)  reduces  to 

Z 

f  (l,0 )  =  -  -  COS  0  =  f (0 ) 

/s 

which  is  no  longer  dependent  upon  displacement  x  and  there¬ 
fore, 

f'(x,0)  =  0. 

The  function  g(x,0)  remains  the  same  as  in  equation  (44). 
Therefore,  equation  (22)  is  reduced  to  the  condensed  form 

wlv(x,0)  +  w(x,0)  +  f(0)  w"(x,0)  =  g(x,0).  (47) 

As  a  check  on  the  validity  of  the  form  of  equation  (47) 
the  uncoupled  equations  of  motion  ( 37 *  )  and  (38* )  can  once 
again  be  derived  from  equation  (47)  by  using  the  technique 
of  normal  mode  expansion. 

If  the  same  finite-difference  scheme  is  psed  on  equa¬ 
tion  (47)  as  was  used  on  equation  (44),  then  the  form  of  the 
recursion  matrix  formula  (46)  remains  the  same  except  that 
the  5^  matrix  of  coefficients  is  now  only  time  dependent 
and  no  longer  depends  upon  displacement,  except  for  dimen¬ 
sion.  This  matrix  becomes  the  time  dependent,  symmetric, 

(n  x  n)  matrix 
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Ir  \r 

vhere  y  and  6  are  the  same  as  those  previously  defined 
except  that  fk  replaces  f  ^ 

—Jr 

With  the  substitution  of  the  above  matrix  for  B  into 
equation  (46),  this  recursion  formula  can  be  programed  into 
the  computer  to  obtain  the  solution  to  the  difference  equa¬ 
tion  representing  the  equation  of  motion  (22). 
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SUMMARY  AND  RESULTS 

The  difference  equations  representing  the  fourth  order 
form  (20)  and  (22)  of  the  equation  of  motion,  and  the 
uncoupled  forms  (37'),  (38')  and  (42)  of  the  equation  of 
motion  were  programed  and  solved  on  a  high  speed  digital 
computer  using  numerical  techniques.  Discounting  diagnostic 
and  de-bugging  time,  the  solution  of  either  difference 
equation  required  approximately  0.87  hours  of  computer  time 
and  the  solution  of  the  first  uncoupled  equation  required 
approximately  0.05  hours  of  computer  time.  In  both  cases 
the  selected  Increment  of  angular  crank  input  was  0.5°  and 
for  the  solution  of  the  difference  equation  50  displace* 
ment  Increments  were  used. 

The  effect  of  varying  the  size  or  number  of  the  A  0  or 
A*  increments  on  a  solution  was  investigated.  It  was 
determined  that,  in  the  solution  of  the  uncoupled  equations 
a  change  from  £  0  =  0.5*  to  A  8  =  0.25*  had  little  or  no 

effect  on  the  solution  to  eight  decimal  places.  In  the 
solution  of  the  difference  equation,  however,  an  increase 

in  A  6  from  0.5*  to  5*  caused  the  solution  to  become 
unstable  and  no  correlation  was  achieved.  It  was  also 
determined  that,  if  the  total  number  of  displacement 
Intervals  is  reduced  from  50  to  10,  then  the  solution 
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of  the  difference  equation  approaches  the  solution  of 
the  first  uncoupled  equation.  See  Figure  7.  Ho  further 
Investigation  of  this  phenomenon  was  attempted. 

Effect  of  Assuming  Variables  Separable 
To  arrive  at  the  uncoupled  equations  of  motion  from 
the  simplified,  linear  partial  differential  equation  of 
motion  a  technique  based  on  the  method  of  separation  of 
variables  was  used  (see  expression  (26)).  in  previous 
work  In  the  field  of  vibrating  beams  this  technique  has 
been  successfully  used  to  great  advantage;  however,  for 
this  specific  problem  this  method  falls  short  of  giving 
a  complete  and  accurate  description  of  the  solution  of 
the  governing  partial  dlff erential  equation  because  of 
approximations  needed  In  the.  development  of  the  uncoupled 
equations.  To  compare  the  solution  of  the  difference 
equation  with  the  solution  of  the  first  uncoupled  equation 
see  Figure  6.  From  the  Inspection  of  Figure  6  It  can  be 
seen  that  good  correlation  between  the  solutions  is 
achieved  during  the  first  and  third  revolutions  of  the 
crank  but  during  the  second  revolution  the  two  solutions 
are  completely  out  of  phase  and  the  large  negative  dis¬ 
placement  predicted  by  the  uncoupled  equation  Is  not 
predicted  by  the  governing  partial  differential  equation. 

It  can  therefore  be  seen  that  the  uncoupled  equation 
predicts  a  more  vigorous  vibration  than  does  the  difference 


*7 


o;^ 


&<■ 

VS 


o  o  o  y 

C\J  if)  O  ; 

O  —  OJ 

. . 

CON  w  J 


(e‘x)M-3anindwv  noiivugia 

Figure  6.  Solutions  of  First  Uncoupled  Equation  and  Simpli¬ 
fied  Equation  (22)  50  Displacement  Increments. 
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equation.  The  solution  of  the  difference  equation,  however, 
exhibits  more  Irregularities  and  Is  not  quite  as  repetitive 
In  nature  as  Is  the  solution  of  the  first  uncoupled  equation. 
This  would  tend  to  Indicate  that  the  uncoupled  solution 
Is  only  a  periodic  approximation  of  the  difference  equation. 
The  addition  of  more  uncoupled  solutions,  second,  third, 
etc.,  will  have  little  affect  on  alleviating  this  problem 
of  non-correlation. 

This  Indicates  that  this  method  of  normal  mode  expan¬ 
sion  must  be  used  with  caution  for  this  type  of  problem, 
and  that  Its  validity  depends  upon  the  form  of  the  forcing 
function  (right  hand  side  of  equation  (22)).  In  this 
particular  case,  this  technique  is  not  the  best  because 
the  forcing  function  la  a  linear  function  In  displacement 
and  does  not  lend  itself  to  proper  separation.  During  the 
process  of  uncoupling  this  forcing  function  was  reduced 
to  a  function  of  time  only  by  evaluating  the  integral  on 
the  right  hand  side  of  expression  (29).  In  doing  this  the 
effect  of  the  x  varying  function  is  approximated  by  its 
behavior  at  x  =  L/2  only.  Since  the  right  hand  side 
of  the  equation  of  motion  represents  the  transverse 
component  of  rigid  body  acceleration  for  any  point  on 
the  connecting  rod,  then  this  approximation  does  not  seem 
reasonable.  However,  with  a  decrease  in  the  number  of 
displacement  Increments  used  In  the  solution  of  the 
difference  equation,  the  resulting  solution  approaches  that 
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arrived  at  by  solving  the  first  uncoupled  equation.  See 
Figure  7  for  this  comparison.  This  Indicates  that  the  solu¬ 
tion  of  the  first  uncoupled  equation  does,  In  part,  repre¬ 
sent  the  solution  of  the  partial  differential  equation 
and  can,  therefore,  be  used  to  determine  the  gross  effects 
of  varying  S  or  Z. 

From  Figure  8  It  cam  be  seen  that  the  assumption 
^Pn(x)  =  sin 

Is  not  only  valid  but  that  a  further  assumption  Is  In 
order.  This  assumption  Is  that  all  of  the  modes  can  be 
neglected  except  the  first  mode  or,  in  other  words, 
that  the  solution  will  take  the  form 

w(x,t)  =  q(t )  sin  ££  (47) 

rather  than  be  an  infinite  sum.  If  this  expression  Is 
substituted  Into  the  equation  of  motion  and  if  the  right 
hand  side  of  this  equation  Is  approximated  by  the  first 
term  of  Its  half  range  Fourier  sine  series,  then  a  better 
approximation  of  the  solution  of  the  difference  equation 
may  be  obtained.  No  analytical  work  was  done  to  sub¬ 
stantiate  the  above  statement. 

The  validity  of  the  assumption  In  expression  (47) 
can  also  be  substantiated  If  a  comparison  between  the 
natural  frequencies  <i >  n  and  the  Input  crank  speed  w  is 
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Figure  8.  Vibrational  Amplitude  Distribution  Along  the 

Connecting  Rod  From  Solution  Of  Equation  (20). 
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made.  Since  co  Is  the  angular  frequency  of  the  forcing 
function  a  resonant  condition  will  be  established  for 
any  mode  whose  natural  frequency  Is  close  to  u>,  and  little 
or  no  detectable  excitation  will  occur  for  natural 
frequencies  above  this  value.  The  first  natural  frequency 
can  vary  In  value  above  and  below  the  value  of  the 
constant  u>  depending  upon  the  relative  magnitudes  of  Z 
and  ffaS,  but  the  second  natural  frequency  ,  which  Is 
always  much  greater  than  the  maximum  value  of  «i>  i ,  will 
not  be  excited  by  u>  to  as  great  an  extent  as  u  i .  For 
the  third  and  higher  natural  frequencies  this  Is  also 
true  because  these  frequencies  Increase  In  proportion  to 
n*  so  that  the  excitation  of  the  higher  modes  by  w  Is 
negligible.  This  line  of  reasoning  was  also  substantiated 
during  the  solution  of  the  uncoupled  equations  of  motion 
(S  =  0.2,  Z  =  1 .5,  r  =  2.0  inches  and  L  =  12.0  Inches)  when 
it  was  calculated  that  the  maximum  value  of  vibrational 
amplitude  of  the  second  mode  equation  was  always  less  than 
10" ‘times  the  value  of  the  first  mode  amplitude  and 
that  the  maximum  value  of  the  vibrational  amplitude  for 
the  third  mode  was  approximately  10“a  times  that  of  the 
first  mode. 

Effect  of  Assuming  Small  Slops 
The  assumption  of  small  slope  (implies  small  deflec¬ 
tion  also)  Is  quite  common  In  vibrational  analysis. 
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Figure  9  Is  a  graphical  example  of  why  this  Is  a  valid 
assumption.  This  figure  is  a  graphical  representation 
of  the  solution  of  the  finite  difference  approximations 
for  the  equations  of  motion  (20)  and  (22).  The  only 
difference  between  equations  ( 20 )  and  (22)  is  that  an 
additional  term  proportional  to  the  slope  was  assumed  to 
be  zero  in  the  derivation  of  equation  (22). 

In  order  to  determine  analytically  why  this  assumption 
is  valid  first  compare  the  orders  of  magnitude  of  the 
components  in  the  coefficient  of  w”  in  equation  (43).  It 
can  be  seen  that,  for  reasonably  large  ranges  in  value 
of  the  parameters  S  and  Z,  the  term 

7§~  008  0 

is  the  dominant  factor.  The  assumption  made  above  only 
eliminates  the  remaining  negligible  terms.  Another  term 
eliminated  by  this  assumption  is  the  entire  w'  term  in 
equation  (43).  This  elimination,  at  first,  would  seem 
unreasonable  except  that,  in  the  solution  of  a  partial 
differential  equation,  the  higher  order  terms  have 
the  greatest  effect.  Therefore,  since  the  equation  in 
question  is  a  fourth  order  partial  differential  equation, 
the  fourth  order  term  will  have  the  greatest  effect  on 
determining  the  solution  while  the  first  order  term  will 
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hare  the  least  effect.  This,  of  course,  places  magnitude 
restrictions  on  the  coefficient  of  w*  which,  in  this  case, 
must  have  been  met. 

The  above  is  also  pointed  out  in  examining  the  dif¬ 
ference  equation  (45)  thereby  noting  that  the  only  change 
in  the  recursion  formula  (46)  In  going  from  the  difference 
equation  representing  (20)  to  that  representing  (22)  is 
in  the  matrix.  Therefore,  if  « j  <  <  1  and  f^^  Z  fjj., 
which  is  true  in  this  case,  then  the  solutions  should  be 
close  and,  therefore,  the  assumption  valid. 

gffect  of  Damning 

The  uncoupled  equations  of  motion  (37* )  and  (38') 
were  modified  to  include  damping  and  equation  (42)  re¬ 
sulted.  The  value  of  the  damping  ratio  17  was  increased 
in  large  increments  until  the  effects  of  the  initial 
conditions  on  the  solution  were  eliminated.  It  was  de¬ 
termined  that  a  value  of  r)  =  0.10  was  sufficient  to  cause 
this  phenomenon  and  Figure  1 0  resulted  from  the  solution 
of  equation  (42).  Figure  10  is  a  graphical  comparison 
between  a  damped  and  undamped  solution  on  the  first  un¬ 
coupled  equation.  Is  would  be  expected,  this  figure 
predicts  that  the  mechanism  should  behave  more  sluggishly 
than  is  predicted  when  the  natural  damping  is  assumed  to  be 
negligible.  It  can  be  seen  that  the  peak  amplitudes  occur 
at  approximately  the  same  time  in  the  cycle  but  are  some- 
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what  reduced  in  magnitude.  See  Table  1.  This  fact  may 
indicate  that  the  value  v  =  0.10  is  too  high  an  estimate 
for  the  damping  ratio.  The  important  fact  to  note  is  that 
the  smaller  vibrations  between  the  larger  positive  and 
negative  amplitude  peaks  are  rapidly  damped  out  whereas 
the  large  amplitude  peaks  are  only  partially  affected. 

^Ajcimum  Variations  in  Piston  Displacement 

Table  1  presents  the  maximum  variations  in  piston 
displacement  at  the  four  predominant  vibrational  peaks 
that  occur  during  the  first  two  revolutions  of  the  crank. 

The  amplitude ,  w,  and  the  crank  angle,  0,  were  taken  from 
the  data  represented  in  Figures  6,  9  and  10,  and  Appendix 
C  contains  the  derivation  for  the  expressions  representing 
A  L  and  A  X0. 

It  can  be  noted  that,  in  the  case  of  Table  1 ,  the 
deflection,  w,  is  considerably  larger  than  the  piston 
variation,  A  Xfi.  This  is  due  to  the  fact  that  A  1  is 
proportional  to  the  square  of  the  amplitude,  a  small  number, 
and  inversely  proportional  to  the  length  of  the  connecting 
rod. 

For  the  first  two  amplitude  peaks  all  of  the  solutions 
itemized  in  Table  1  predict  comparable  variations  at  the 
same  time  in  the  cycle.  During  the  second  cycle  divergence 
in  both  piston  variation  and  time  of  occurrence  beoomes 
pronounoed • 
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Table  1 .  Predicted  Variation  In  Piston  Displacement 


(r  =  2 

in,  L 

=  12  in,  S  = 

0.2,  Z  =  1 . 

5) 

Main 

Peak 

Equation 

0 

w 

Amplitude 

Inches 

A  L 

Shortening 

Inches 

Inches 

1 

(20) 

155* 

.0904 

.00168 

.00167 

(22) 

160* 

.0981 

.00198 

.001968 

(37* ) 

135* 

.0908 

.001698 

.001688 

(42) 

150° 

.0889 

.001625 

.001615 

2 

(20) 

250° 

-.167 

.00573 

.00565 

(22 

255° 

-.180 

.00666 

.00657 

(37  ) 

245° 

-.157 

.00507 

.005 

(42) 

240° 

-.141 

.00408 

.00402 

3 

(20) 

455° 

.  .158 

.00507 

.00493 

(22) 

460° 

.173 

.00615 

.00597 

(37  ) 

505° 

.178 

.00651 

.00646 

(42) 

510° 

.0889 

.001625 

.001615 

4 

(20) 

630° 

-.125 

.00321 

.00312 

C22) 

635® 

-.1097 

.002475 

.0024 

(37*) 

595° 

-.207 

.0088 

.00867 

(42) 

600* 

-.141 

.00408 

.00402 

Discussion  of  the  Parameters 
In  order  to  discuss  the  possible  effeot  that  the 
dimensionless  parameters  S  and  Z  will  have  on  the  vibra¬ 
tional  characteristics  of  a  mechanism  a  sample  slider-crank 
mechanism  will  be  analyzed. 

Equations  (39)  and  (40 )  are  expressions  representing 
S  and  Z  respectively.  If  a  value  of  Z  =  1.0  is  chosen  for 
the  physical  parameter  and  if  the  linkage  properties  are 
reasonably  chosen  then  the  following  quantities  can  be 


59 


specified 

I  =  30  x  10*  lb/ln2,  r  =  1  in, 

I  =  0.0208  in4  (1  x  i  L  =  6  in, 

cross-section), 

M3  =  0.0022  lb-sec*/in,  H*  =  0.011  lb-seca/ln. 

The  S,  or  speed,  parameter  cannot  be  specified  at  this 

time  because  the  operating  input  speed  u>  of  the  mechanism 
has  not  been  specified.  These  quantities  will  be  used 
throughout  this  subsection  for  illustration  purposes. 

In  the  uncoupled  equations  of  motion  (371 )  or  (38' ) 
the  coefficient  of  qn(0)  represents  the  natural  frequencies 
of  the  modes  of  vibration.  Written  in  terms  of  the 
dimensionless  parameters  this  coefficient  becomes 


[ff*n4S  +  ff*nBZ  cos  8) 


(48) 


Since  0  ‘6  s  ",  the  second  term  in  this  expression  will 
become  negative  when  the  angle  6  is  in  the  second  or 
third  quadrant  and  Cos  6  s  -1.0  for  8  equaling  ir  or  any 
odd  multiple  of  rr. 

It  can  therefore  be  shown  that  the  uncoupled  equations 
of  motion  no  longer  predict  the  vibratory  motion  of  the 
connecting  rod  unless 

Z  *  n*  S  (49) 

for  any  given  situation.  If  Z  exceeds  this  limit  then 
the  first  natural  frequency,  as  calculated  from  expression 


60 


(48),  will  be  imaginary,  a  situation  that  is  contrary  to 
physical  reality. 

It  is  interesting  to  note  that  the  criteria  stated  in 
equation  (49)  corresponds  to  the  situation  where  the  dynamic 
axial  load  approaches  the  Euler  criteria  for  static  buckling 
of  long,  slender  columns.  That  is 

K  Cos  0  =  Pcr 

where 


ira  El 


Although,  in  general,  the  connecting  rod  will  not  have  a 
sufficiently  low  slenderness  ratio  to  allow  the  use  of 
the  Euler  buckling  criteria  directly,  this  condition  is 
a  guide  and  loadings  above  Pcr  should  be  avoided  by  the 
designer. 

S  and  Z  are  the  parameters  that  govern  the  behavior  of 
the  equation  of  motion.  However,  since  the  equation  of  motion 
has  not  been  completely  nondimensionallzed,  two  other  para¬ 
meters  of  the  mechanism  must  be  specified  in  order  to  attempt 
a  solution;  these  are  the  crank  length,  r,  and  the  con¬ 
necting  rod  length,  L.  Since  the  parameters  S  and  Z  are 
Important  in  analyzing  the  vibratory  properties  of  a  slider-crank 
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mechanism,  a  discussion  of  the  effect  of  mechanism 
properties  on  these  parameters  is  necessary. 

First  consider  only  the  dimensionless  ratio  r/L.  If 
r/L  remains  constant  then  Z  will  remain  constant;  however, 
this  can  be  true  even  though  r  and  L  both  vary  in  a  manner 
consistent  with  the  constant  r/L  ratio.  If  r  increases 
then  L  must  increase  proportionately.  This  increase  in 
L  will  cause  a  large  decrease  in  S,  since  S  is  inversely 
proportional  to  L3 •  This  has  the  effect  of  making  s*S 
approach  Z,  and,  as  will  be  shown  in  detail  later,  as  these 
two  quantities  become  equal  in  value  the  vibrational 
oscillations  of  the  connecting  rod  become  greater  in  both 
number  and  amplitude.  This  has  a  very  decided  effect  if 
Z  is  .close  to  ir8S  before  L  is  increased  because  a  slight 
change  in  S  in  this  range  has  a  very  marked  influence  on 
the  amplitude  of  vibration.  The  corresponding  increase 
in  r  also  causes  an  increase  in  the  non-homogeneous  part 
of  the  equation  of  motion;  however,  the  effect  of  in¬ 
creasing  L  is  much  more  pronounced  than  is  the  proportionate 
increase  in  r.  The  converse  of  this  statement  is  also 
true. 

Next  consider  the  mass  ratio  H,/Ma  .  If  this  ratio 
remains  constant  then  Z  will  remain  constant,  but 
both  M*  and  M,  can  vary  as  long  as  the  specified  ratio 
is  maintained.  If  Ms  increases  then  S  decreases,  since 
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S  is  inversely  proportional  to  M,  ,  and  tr*S  approaches  Z. 
This  condition  will  cause  an  increase  in  vibrational  ac¬ 
tivity.  However,  since  ^  has  less  of  an  effect  on  S  than 
does  1,  the  vibrational  effect  mentioned  in  the  above 
paragraph  will  be  much  more  pronounced. 

The  modulus  of  elasticity,  E,  depends  only  upon  the 
material  used  in  the  connecting  rod,  and  the  moment  of 
inertia,  I,  depends  upon  the  cross-sectional  area  of  the 
connecting  rod.  il though  the  freedom  with  which  I  can  be 
varied  is  limited  and  the  variation  of  I  is  somewhat 
limited  by  the  mass  Ma ,  these  two  parameters  can  be  varied 
independently  of  Z,  Increasing  the  El  combination  in¬ 
creases  S  and,  therefore,  decreases  the  vibrational  oscilla¬ 
tions.  The  converse  of  this  statement  Is  also  true. 

The  only  remaining  quantity  that  has  an  effect  on 
the  relative  value  of  the  dimensionless  parameters  S  and 
Z  is  the  angular  input  speed  of  the  crank,  go.  The  speed 
parameter  S  (This  parameter  was  given  this  name  because 
it  is  an  effective  measure  of  the  angular  velocity  of  the 
crank.)  is  inversely  proportional  to  go  a.  It  must  be 
noted  that  go  is  the  only  term  in  S  that  can  be  varied  over 
a  wide  range  of  values  without  affecting  the  remainder 
of  the  mechanism  parameters. 

Is  would  be  expected  when  go  =  0  the  speed  parameter 
S  =  ®.  This  in  turn  predicts,  through  the  equation  of 
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motion,  that  no  vibration  is  present  in  the  connecting 
rod.  Is  S  decreases  (or  «  increases)  the  inequality 
stated  in  expression  (49)  approaches  the  limiting  oase 
Z  ss  saS 

The  value  of  w  that  gives  this  equality,  for  a  speci¬ 
fied  mechanism,  is  the  predicted  maximum. value  for  the 
input  speed  of  the  crank.  As  an  example,  using  the 
quantities  given  earlier  . 


%ax  = 


ffaBI 

H,LS 


=  3600  rad/sec 


=  34,300  RPM. 


To  further  illustrate  the  effect  of  increasing  r 
and  L  while  keeping  the  ratio  r/L  constant,  let  r  »  2  in. 
and  L  =  12  in.  Performing  the  calculations  for  u^ax 
yields 

c^ax  =  1273  rad/sec  =  12,150  RPM. 

It  can  be  seen  that  the  maximum  allowable  angular  input 
is  greatly  reduced  by  this  increase  in  L. 

However,  if  Ha  is  doubled  while  the  mass  ratio  is 
held  constant  (The  same  would  be  true  if  E  or  I  would  be 
decreased  by  one  half.)  then  this  calculation  yields 

Snax  =  2550  rad/sec  =  24»400  HPM* 

The  effect  of  this  change  on  the  maximum  permissible 
running  speed  is  much  less  than  that  recorded  from  the 
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change  In  connecting  rod  length* 

The  effects  of  increasing  the  angular  speed  co  on  the 
predicted  vibrational  characteristics  of  a  mechanism  are 
recorded  in  Figure  1 1 .  The  data  for  the  figure  were 
generated  by  solving  the  first  uncoupled  equation  of  motion 
with  r  =  2  in.,  L  =  1 2  in.,  Z  =  1 .0  and  0.13  *  3  «  1.0. 
Using  the  necessary  quantities  from  those  given  at  the 
beginning  of  this  subsection  the  parameters  can  be  trans¬ 
lated  into  representing  the  following  variations  in  speed 


W 


» 


therefore 

w  =  405  rad/sec  ~  S  =  1 .0, 

574  rad/sec  ~  0.5» 

906  rad/sec  ~  0.2, 

1047  rad/sec  ~  0.1 5« 

1123  rad/sec  ~  0.13* 

It  can  therefore  be  seen  that  for  speeds  that  are  slow  in 
comparison  to  the  maximum ^predic ted  speed,  small,  smooth 
vibrations  are  recorded  with  amplitude  reversals  occurring 
only  every  180*.  is  the  speed  is  increased  the  amplitude 
of  vibration  increases  but  the  reversals  still  occur  only 
every  180°.  Increasing  the  speed  further,  0.15  *  S  *  0.2, 
finally  causes  the  reversals  to  become  more  numerous.  At 
this  speed,  and  above,  the  vibrations  increase  in  both 
amplitude  and  number  at  an  accelerated  rate  until  the 
equation  of  motion  becomes  unstable  and  large  unbounded 
vibrations  are  predicted. 


Z=I.O 
r-2.0  IN. 


CD 

i 
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Figure  11.  Solution  of  First  Uncoupled  Equation  for 
Different  Crank  Speeds. 
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Although  not  recorded  on  Pigure  1 1  the  first  uncoupled 
equation  was  also  solved  for  S  =  0.11,  a  value  that  barely 
exceeds  the  minimum  permissible  value  of  S.  The  trend 
recorded  In  Figure  1 1  was  continued  and  a  peak  amplitude 
of  approximately  0.9  was  reoorded  at  0  =  560® ;  this 
represents  a  piston  variation  of  approximately  0.166  Inches. 
This  solution  was  not  Included  In  Figure  1 1  because  unstable 
predictions  were  recorded  during  the  third  revolution  of 


the  crank 
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CONCLUSIONS 

Equations  of  motion  for  the  transverse  vibrational 
characteristics  of  a  slider-crank  mechanism  have  been  derived, 
and  solutions  of  three  different  forms  of  these  equations 
have  been  obtained  for  a  mechanism  with  a  crank  length  of 
2  inches.  The  data  from  these  solutions  are  presented  in 
Figures  6,  8,  and  9.  The  primary  objectives  set  forth  at 
the  beginning  of  this  study  have  been  achieved.  However, 
due  to  the  required  extended  use  of  computer  facilities,  the 
effect  of  extensive  mechanism  parameter  variation  was  not 
obtained. 

The  uncoupling  procedure  used  during  this  investigation 
is  not  entirely  adequate  in  that  it  only  results  in  a  peri¬ 
odic  approximation  of  the  solution  to  the  more  complex  fourth 
order  form  of  the  equation  of  motion.  However,  since  the 
uncoupled  equations  are  sufficient  to  predict  the  gross  vi¬ 
brational  results,  and  since  these  equations  predict  vibra¬ 
tional  activity  greater  than  that  predicted  by  the  linear 
fourth  order  form  of  the  equation  of  motion,  the  uncoupled 
equations  can  be  used  for  design  purposes.  In  light  of  the 
results  recorded  in  Figure  8  a  better  approximation  of  the 
solution  may  be  obtained  by  assuming  that  the  amplitudes  of 
all  of  the  modes  of  vibration  are  negligible  except  the  first 
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mode.  The  mode  represented  hy  the  first  natural  frequency 
is  the  predominant  mode  because  the  value  of  the  frequency 
of  the  forcing  function  is  close  to  the  value  of  the  first 
natural  frequency.  The  second  and  higher  natural  frequencies 
are  always  much  greater  than  the  first  natural  frequency, 
therefore,  excitation  of  these  modes  is  much  less. 

If  the  restrictions  placed  on  the  analysis  by  the  load¬ 
ing  approximations  derived  in  Appendix  B  are  too  restrictive, 
then  the  solution  of  equation  ( 1 6 * )  can  be  undertaken.  Ob¬ 
taining  this  solution  would  follow  the  same  procedure  as 
obtaining  the  solution  of  equation  (20 ).  The  only  differ¬ 
ence  will  occur  in  the  complexity  of  the  5^  matrix  in  the 
recursion  formula. 

Figure  9  shows  excellent  correlation  between  a  more 
exact  and  a  simplified  version  of  the  linear  fourth  order 
equation  of  motion.  Appropriate  simplifying  assumptions  can, 
therefore,  be  made  without  destroying  the  validity  of  the 
resulting  solution.  However,  assumptions  should,  if  pos¬ 
sible,  be  investigated  to  assure  their  appropriateness. 

As  pointed  out  in  Figure  10  a  relatively  small  amount 
of  damping  does  not  affect  the  major  vibrational  amplitude 
peaks  to  any  great  degree.  Damping  does,  however,  markedly 
affect  the  number  of  complete  amplitude  reversals  so  that 
minor  vibrational  amplitude  peaks  are  rapidly  damped  out. 

In  evaluating  the  effects  of  parameters  on  the  vibra¬ 
tional  characteristics  of  the  mechanism  it  is  discovered 
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that  the  lengths  r  and  L,  especially  L,  have  a  far  greater 
effect  on  the  response  of  the  system  than  do  other  physical 
parameters  such  as  type  of  material,  cross-section  or  mass. 
Increasing  L  causes  a  marked  Increase  In  vibrational  activity 
and  also  greatly  reduces  the  maximum  operating  Input  speed. 

To  completely  evaluate  the  effect  of  parameter  varia¬ 
tion  on  the  transverse  vibration  of  the  connecting  rod, 
equations  (37'),  (38*)  and  (43)  should  be  non-dimensionallzed 
as  suggested  and  solved  for  various  values  of  S,  Z  and  r/L. 

Prom  this  Investigation  the  designer  is  give  one  more 
tool  with  which  to  work.  Noting  that  the  vibrational  ampli¬ 
tude  increases  as  7raS  approaches  Z,  if  accuracy  of  piston 
displacement  is  required  then  the  designer  must  avoid  the 
situation  where  Z  is  approximately  equal  to  it9- S.  To  do  this 
he  has  at  his  disposal  both  the  variation  of  any  physical 
linkage  parameter  and  the  variation  of  the  crank  speed. 

Based  on  the  results  in  Figure  11,  if  iraS  is  approximately 
twice  the  value  of  Z  then  the  predicted  transverse  vibrations 
rapidly  become  negligible. 

Certain  assumptions  were  made  during  this  vibrational 
investigation  to  arrive  at  an  equation  of  motion.  The  effect 
of  the  majority  of  these  assumptions  cannot  be  adequately 
evaluated  analytically .  Therefore,  this  analysis  should  be 
supplemented  with  an  experimental  effort. 
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APPENDIX  A 

Rigid  Body  Acceleration  Analysis  of  a  Slider-Crank  Mechanism 
In  Figure  1  the  length  of  the  line  segment  X75  =  r  sin  0; 
the  length  of  line  segment  =  r  eos  0,  and  the  length  of  the 
line  segment  755  =  (La  -  ra  slna  0)^.  The  well  known  equa¬ 
tion  relating  the  piston  displacement  to  the  angular  Input 
lsf  therefore,  derived  by  noting  that 

X0  =  57  +  755  =  r  cos  0  +  L(  1  -  ra/La  slna  0)^.  (1) 

By  making  use  of  the  binomial  expansion  and  by  restricting 
the  Investigation  to  small  r/L  ratios  (r/L  <  £),  the  radical 
term  can  be  replaced  by 

(1  -  ra/La  slna9)^  =  1  -  i  ra/La  slna  0. 

Therefore, 

XB  =  r  cos  0  +  L(1  -  ra/La  slna  0).  (la) 

The  approximate  acceleration  of  the  piston,  B,  Is  found 
by  differentiating  equation  (la)  twice  with  respect  to  time. 

XB  =  -  r  0  sin  0  -  r*/2L  0  sin  20, 

XB  =  -(r  cos  0  +  ra/L  cos  20)0®  -  (r  sin  8  +  ra/2L  sin  20)0. 

(2a) 
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The  x  and  y  components  of  the  rigid  body  acceleration 
of  a  general  point  4  on  the  connecting  rod  are  determined  by 
summing  the  appropriate  components  of  the  acceleration  of 
point  A  and  of  the  relative  acceleration  of  4  with  respect 
to  A, 


XA  = 

r  cos 

9, 

YA  = 

r 

sin  0, 

xa  = 

-  r  0 

sin 

9, 

e 

ya  = 

r 

• 

0  cos  0, 

•  • 

XA  = 

-  r  0 

cos 

0  - 

•  •  *  • 

r  0  sin  0,  YA  = 

- 

e  •• 

r  0  sin  9  +  r  0  cos  0, 

Using  an  orthogonal  coordinate  transformation  from  the  (X,Y) 
system  to  the  (xfy)  system  gives  the  components  of  the  accel¬ 
eration  of  point  A  In  the  x  and  y  directions 

• «  •  •  *  * 

xA  =  XA  cos  0  -  Ya  3ln  <o. 

«  •  a  e  •  * 

yA  =  XA  sin  to  -+  Ya  cos  <n. 

The  acceleration  of  some  arbitrary  point  4  on  the  connecting 
rod  can  now  be  written  as 
«U.t)  =  =  xk  +  A€a 

•  *  e  • 

a  XA  C08  <0  -  Ya  sin  (0  -  4  <0a 

=  -  r0a(cos  0  cos  <p  -  sin0  sin  <p) 

•  •  • 

-r0(sln  0  cos<0  +  cos  0  sin  <P)-  4  <Pa ,  (3a) 

=  V  |  =  VA  +  a|a 

=  XA  sin  +  Ya  cos  V  +  i 

e 

=  -  r0(cos  0  sin  <o  +  sin  0  cos  0) 


-  r0(sln  0  sin  -  cos  0  cos  <P)  +  4 


(4a) 
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In  equation  (4a)  and  (3a)  A3^  represents  the  normal  compon¬ 
ent  of  the  acceleration  of  Z  relative  to  A  and  A^a  represents 
the  tangential  component  of  the  acceleration  of  Z  relative 

•  e 

to  A.  (Positive  9  must  be  counterclockwise  about  the  point 
A  for  this  to  hold.) 

Prom  the  coordinate  transformation  it  can  also  be  shown 


that 

••  .#  • 

9(L,t)  =  Xjj  sin  to 

•  e  •• 

a(L,t)  =  XB  cos  9. 

The  above  equations  have  served  as  a  check  on  the  derivation. 

Two  relationships  that  can  be  used  to  express  9  as  a 
function  of  9  are  also  available  from  the  geometry  of  Figure 


1 


sin  <b  =  r/L  sin  9, 

cos  <e  =  (1  -  ra/La  sina  9)^. 


(5a) 

(6a) 


e  e  e 

Bquation  (5a)  can  be  used  to  determine  to  and  <o  as  functions 
of  the  angular  input  9.  The  first  derivative  of  (5a)  with 
respect  to  time  yields 


(7a) 


Care  must  be  taken  to  assure  that  the  sign  of  this  function 
representing  9  agrees  with  the  assumed  convention  of  counter¬ 
clockwise  rotation  about  A.  It  can  be  seen  from  a  velocity 
polygon  representing  the  velocities  of  the  mechanism  in  the 


74 


position  shown  in  Figure  1  that  the  relative  velocity  of  B 
with  respect  to  A  will  be  directed  in  the  negative  y  dlrec- 
tion.  Therefore,  to  keep  <p  within  the  chosen  algn  convention 


(8a) 


Differentiating  equation  (7a)  with  respect  to  time  yields 
the  angular  acceleration  of  the  connecting  rod 

•  '  *  «*  0  •* 

-sin  <p  <pa  +  cos  0  <p  =  £(-  sin  0  8'  +  cos  0  0), 

*  =  I' cos'^(“8ln9  9*  +  008  9  9)  +  cosHp  **•  (9a) 

Substituting  equation  (8a)  into  equation  (9a)  gives  the  ex¬ 
pression 


e  • 

= 


cos  <p 


(-sin  0  0a  +  cos  00)+  (£^) 


a  casa0slm> 
COS3V>  « 


It  can  be  seen  from  an  acceleration  polygon  representing  the 
accelerations  of  the  mechanism  in  the  position  shown  in  Fig¬ 
ure  1  that  the  tangential  component  of  the  acceleration  of  B 
relative  to  A  is  directed  in  the  negative  y  direction  which 
means  that  ‘<p  is  clockwise  about  A.  Therefore,  the  keep*<p 
within  the  chosen  convention  of  positive  being  counterclock¬ 
wise  a  negative  sign  must  be  included.  Therefore, 

'i  -  rdfr?  (sln  6  •*  - 008  9  •>  -  <?>■ 

If  the  angular  input  0(t)  is  a  linear  function  of  time, 
0  =  w  t,  then  the  basic  equations  given  above  reduce  to 
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Xg  =  -  r  toa(cos  0  +  £  cos  20), 

a£4,t)  =  -  r  toa(cos  0  cos  <p  -  sin  0  fcin  <p)-4  <P* , 

•  • 

/3U,t)  =  -  r  co a( cos  <j  sin  q  +  sin  0  cos  <p)*i  <t>. 


(2a’) 

(3a* ) 
(4a*) 


(8a') 

(10a') 


It  oan  readily  be  seen  that  a  and  /3  are  linear  functions 
in  displacement  and  for  simplicity  can  be  written  as 


aU,t)  =  gi(t)  +  igait)  (5a  ) 

V(  4 ,  t )  =  ga(t)  +  Cg4(t) 

where  the  functions  of  time,  or  0,  are  described  in  equations 
(3a' ),  (4a'),  (8a' )  and  (10a')  above. 
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APPENDIX  B 

Rigid  Body  Force  Analysis 

If  the  assumption  Is  made  that  the  vibration  of  the 
connecting  rod  does  not  affect  the  pin  reactions  at  A  and  Bf 
then  Rx  and  By.  become  known  functions  of  time,  8(t),  that  can 
be  determined  directly  from  equations  (13)  and  (11)  respec¬ 
tively.  These  functions,  however,  are  relatively  complex 
In  nature  and  do  not  lend  themselves  to  simple  mathematical 
formulation.  With  the  hope  of  determining  a  simple  approx¬ 
imation  for  Rx  and  Ry  that  can  be  used  In  equation  (14)  there¬ 
by  eliminating  the  non-linearities  of  the  final  equation  (16), 
a  FORTRAN  program  for  the  IBM  7090  was  compiled.  This  pro¬ 
gram  was  used  to  calculate  the  components  of  the  A  and  B  pin 
reactions  perpendicular  and  parallel  to  the  connecting  rod. 

The  slider-crank  mechanism  analyzed  was  assumed  to  have  a 
constant  Input  crank  speed,  w;  the  connecting  rod  was  assumed 
to  be  a  long,  thin  rod,  and  the  external  piston  force,  P(t), 
was  assumed  to  be  zero.  The  dimensionless  ratios  r/L  and 
Ma/M«  were  varied  to  determine  the  effect  that  these  ratios 
had  on  the  reactions. 

The  non-dlmenslonallzed  results  of  this  program  are 
graphically  illustrated  in  Figures  12,  13,  14  and  15»  The 
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Figure  12.  Pin  Reaction  Rx  Compared  with  -  K  cos  0 
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Ry  Compared  with  -  Q  sin  0 
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-3.0 


approximations 


(17) 


=  -  K  cos  0, 

Ry.  =  -  Q  sin  0, 

seem  to  be  valid  when  the  r/L  ratio  is  small  (r/L  <  J)  and 
the  Ms/M*  ratio  is  small  (Ms/H*  <  1.5).  The  value  of  the 
positive  oonstant  K  can  be  approximated  by 

K  =  (Ms  +  M*  )  r  uP.  (18) 

The  value  of  the  positive  constant  Q  was  not  required  for  the 
solution  of  the  equation  of  motion  but  can  be  approximated 
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APPENDIX  0 


Derivation  of  Equations  to  Determine 
Variation  In  Piston  Displacement. 


If  the  aseumptlon  Is  made  that  the  length  of  the  beam 
L  remains  constant  then 
L*  L' 

L  =  J  o  d  s  =  f0  +  d  x* 


Substituting  the  binomial  expansion  for  the  radical  and 
neglecting  all  but  the  first  term  yields 


L' 

L  =  J  [1+*(v')a 


+  j-(w'  )4+  - — ]dx 


V 


L*  +  L  *<"')•  d*- 


Since  4  L  :  L  -  L' 

and  since  A  L  will  be  small «  it  can  be  approximated  by 
L 

A  L  =  J  £(w'  )*  dx. 
o 

Substituting  the  derivative  of  expression  (47)  into  the  above 
yields 

L  qir  rrx  naq* 

A  L  =  j  £(-£)*  cos*  -j  dx  =  (1c) 

0 

where  q  is  the  amplitude  w(jf,d)  taken  from  Figures  6,  7,  9» 


10  or  11 
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The  position  of  the  piston  during  a  vibration  can  be 
written  as 

Xg  =  r  cos  0  +  I'  (1  -  (£,  )a  slna 

Since  the  vibration  has  been  assumed  small  r/L  =  r/L* . 
Therefore,  substituting  r/L  into  the  above  equation  and  sub¬ 
tracting  it  from  equation  (1)  yields 

A  XB  =  XB  -  X£  =  A  L(1  -  (£)a  sin*  0)^ 

the  vexiation  in  piston  displacement. 


